Higher moment singularities explored by the net proton non-statistical fluctuations 
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We use the non-statistical fluctuation instead of the full one to explore the higher moment sin- 
gularities of net proton event distributions in the relativistic Au+Au collisions at ^Jsnn from 11.5 
to 200 GeV calculated by the parton and hadron cascade model PACIAE. The PACIAE results 
of mean (M), variance (a 2 ), skewness (S), and kurtosis («) are consistent with the corresponding 
STAR data. Non-statistical moments are calculated as the difference between the moments derived 
from real events and the ones from mixed events, which are constructed by combining particles 
randomly selected from different real events. An evidence of singularity at -^/sjvjv ~ 60 GeV is first 
seen in the energy dependent non-statistical S and Sa. 

PACS numbers: 25.75.Dw, 24.60.Ky, 24.10.Lx 



o 

=5 



> 

m 

in 

in 
o 

(N 



I. INTRODUCTION 

Describing the QCD phase diagram as a function of 
temperature T and baryon chemical potential /is is one 
of the fundamental goals of heavy-ion collision experi- 
ments pj. The finite temperature lattice QCD calcu- 
lation at ^tB=0 predicts that a crossover transition from 
the hadronic phase to Quark Gluon Plasma (QGP) phase 
may occur at temperature of 170-190 MeV How- 
ever, a QCD based model calculation indicates that the 
transition could be first order at large /ib [3 • Once the lo- 
cation of the QCD critical point (QCP), where the phase 
transition proceeds from first order to crossover, is iden- 
tified, the global structure of the phase diagram is known 

MM- 

A characteristic feature of QCP is the divergence of the 
correlation length £ 0] and the extremely large critical 
fluctuations Q. In a static and infinite medium, various 
moments of conserved quantities such as the net-baryon, 
net-charge, and net-strangeness are related to the corre- 
lation length £ [9j. Typically the variance (a 2 ) of these 
distributions is related to £ as a 2 ~ £ 2 Q. It is pointed 
out in [loj that higher moments of conserved quantity 
event distributions, measuring deviations from a Gaus- 
sian, are sensitive to the QCP fluctuations. 

In the first quotation of [ll| the hadronic and quark- 
gluon coexisting phases in a small volume have been in- 
vestigated with a simple effective model. Assuming this 
finite system is in a heat reservoir, the order parameter 
(energy density e) fluctuation near the phase transition is 
studied in a way similar to the Landau theory. A positive 
skewness (S) was predicted and a negative kurtosis (k) 



was expected for the extensive thermodynamical quanti- 
ties. In the second quotation of [ll[ the dynamical change 
of skewness and kurtosis was analyzed during hadroniza- 
tion of QGP. It was shown that the skewness changes 
from negative to positive during the transition. While 
the kurtosis is positive in the initial dominantly QGP 
phase and after hadronization, but it is negative during 
the process of transition. Similarly, it was predicted in 
the effective theory [l(], EH that a crossing of the phase 
boundary may result in a change of sign of skewness as 
a function of the energy density. It was also reported re- 
cently that the sign of kurtosis could be negative as well, 
if the QCP is approached from the crossover side of the 
QCD phase transition [13l |. 

The products of higher moments, such as Sa and na 2 , 
are related to the ratio of conserved quantity number 
susceptibilities (x): Sa ~ x^/x^ an d K,a 2 ~ x^/x^ 
[i"4j . It is predicted that the conserved quantity event 
distribution becomes non-Gaussian and the susceptibility 
diverges when QCP is approached. This causes Sa and 
na 2 to change significantly. 

Recently, QCP and the higher moments of conserved 
quantity event distribution in heavy ion collisions at BNL 
Relativistic Heavy Ion Collider (RHIC) energies have 
aroused further interest both experimentally 0, UH, [l6[ 
and theoretically pJHSlj]- On the experimental side, a 
method to determine T c based on data has been pro- 
posed in [Pa ]. By comparing the lattice results with the 
RHIC BES (Beam Energy Scan) fluctuation data of vari- 
ance, skewness, and kurtosis systematically, the critical 
temperature is determined to be T c = 175*7 MeV. How- 
ever, no evidence of singularity in energy dependence (at 
given centrality) and/or centrality dependence (at given 



2 



0-5 % most central Au+Au @/ s NN = 39 GeV 

0.4<p T <0.8 GeV/c, lyl<0.5 

• STAR 
O PACIAE 



O • 



3 m" r 
S 

ZlO^n 



-°o 

• O 

• o 
• o 

. o 



10 
Net-proton 



FIG. 1: Net-proton event distributions in Au+Au collisions 
at t/snn=39 GeV. The solid and open symbols are STAR 
data (taken from first quotation of Ref. [lg) and the results 
of PACIAE real events, respectively. 



energy) has been reported. On the other hand, copious 
models have been further proposed: such as the Lattice 
QCD [17|, 2+1 flavor Quark- Meson model and Polyakov- 
Quark-Meson mode l fl8| , Polyakov-Nambu- Jona-Lasinio 
(PNJL) model I19L |2CI. Dvson-Schwinger equation |2l| . 

statistical modelpJBl' and the Ur Q MD and AMPT 
transport models [la . l24j etc. Each model has its merits 
and results, but the consistency is lacking and the contra- 
diction is existing among the theoretical models. The dif- 
ficulty is that the fluid dynamical development of energy, 
momentum, and baryon charge density leads to a com- 
plex spatial distribution of the measured and statistically 
analyzed quantities at the freeze-out and hadronization 
domain of the space-time. These also influence the mea- 
sured higher moments even without a phase transition, 
and to disentangle the two effects is not easy 25] . Thus 
the question is still open and further studies are required. 

In this paper a new method is studied to provide more 
insight. Namely, non-statistical moments of net proton 
event distribution are calculated as the difference be- 
tween the moments derived from real events generated 
by the parton and hadron cascade model PACIAE (2(| 
and the ones from mixed events which are randomly con- 
structed according to the real events [13] ■ In this way the 
single particle distribution arising from the fluid dynamic 
development is separated from the two and more particle 
correlations stemming from the PACIAE model where 
the hadronization and freeze-out takes place, according 
to our expectation. 



II. MODELS 



collisions according to the collision geometry and total 
NN cross section. The collision time is calculated for 
each NN collision pair assuming straight line trajectory 
between two consecutive NN collisions and the NN col- 
lision list is then constructed by all collision pairs. A 
NN collision with earliest collision time is selected from 
the collision list and performed by PYTHIA with string 
fragmentation switches-off and diquarks (anti-diquarks) 
break into quark pairs (anti-quark pairs). Thus a parton 
initial state (quarks, antiquarks, and gluons) is eventu- 
ally obtained when NN collision pairs are exhausted. 

The parton rescattering is then proceeded by Monte 
Carlo method using 2 — > 2 leading order perturbative 
QCD cross sections f29|. The parton evolution stage is 
followed by the hadronization at the moment of partonic 
freeze-out (exhausting the partonic collisions) . The Lund 
string fragmentation model and a phenomenological co- 
alescence model are provided for hadronization. After 
this the rescattering among produced hadrons is dealt 
with the usual two body collision model 26]. Only the 
rescatterings among ir, K, p, n, p(uj), A, A, S, S, Q, J/ty 
and their antiparticles are considered for simplicity. 

Like other transport (cascade) models, such as above 
mentioned UrQMD [H and/or AMPT [p], PACIAE 
does not assume equilibrium. It just simulates dynam- 
ically the whole relativistic heavy ion collision process 
from the initial partonic stage to the hadronic final state 
via the parton evolution (not implemented in UrQMD), 
hadronization, and hadron evolution according to a copi- 
ous dynamical ingredients (assumptions) introduced rea- 
sonably. Therefore it is parallel to the experimental 
nucleus-nucleus collisions. These dynamics correctly de- 
scribes the particle, energy, and entropy etc. develop- 
ments, while intensive thermodynamical quantities are 
not defined in this non-equilibrium regime. Messages 
brought by the produced particles in these transport (cas- 
cade) models are all of dynamical origin. Unlike most of 
the hydrodynamic models where the phase transition is 
described via the assumptions for equation of states, in 
PACIAE (the same in UrQMD and/or AMPT) we do 
not implement a phase transition congenitally. However, 
once there is a phase transition signal in transport (cas- 
cade) model calculations, it must be a result of dynamical 
evolution. Of course, further studies are then required. 

The n th moment about the mean of a conserved quan- 
tity, x, with the event distribution P(x) is expressed as 

M<"> = ((.x - (*))") = J (x - (x)) n P(x)dx, (1) 

where the n*' 1 moment about zero reads as 



The PACIAE model [26[ is a parton and hadron cas- 
cade model, which is based on PYTHIA [28]. PACIAE 
consists of the parton initiation, parton evolution (rescat- 
tering) , hadronization, and the hadron evolution (rescat- 
tering) four stages. 

In the parton initiation stage, a nucleus-nucleus col- 
lision is decomposed into binary nucleon-nucleon (NN) 



(x^) = / x n P(x)dx, 



(2) 



The higher moments as well as their products investi- 
gated widely are then 
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FIG. 2: Centrality dependent moments of net-proton event 
distribution in Au+Au collisions at ^/snn~S9 GeV. The solid 
and open circles are STAR data (taken from first quotation of 
Ref. [ly]) and the results of PACIAE real events, respectively. 



skewness: S = M (3) /(Af (2) ) 3/2 , (4) 

kurtosis: k = M^/(M^) 2 - 3, (5) 

5(t = M (3 Vm (2) , (6) 
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FIG. 3: (Color online) Energy dependence of na 2 of net- 
proton event distribution in Au+Au collisions. The solid and 
open circles are STAR data (taken from pj) and the results 
of PACIAE real events, respectively. 



ones from mixed events. Thus the dynamical fluctuations 
are pronounced from the underlying dominant statistical 
fluctuations and are expected to be seen easily in the 
non-statistical higher moments. 

If the n th moment about the mean calculated by real 
(mixed) events is denoted by (M^ ) , then the cor- 

responding n th non-statistical moment about the mean, 
like in 127], is: 



f (n) 
N ON 



r(«) 
L M ■ 



(8) 



and 

no 2 = M^/M^ -3Af( 2 >, (7) 

besides the mean M = and M^> = 0. 

The study of higher moment singularities is a mat- 
ter of dynamics, but the statistical fluctuation is always 
dominant. Therefore, we study the non-statistical part 
of fluctuations instead of full fluctuations in this pap er. 
The parton and hadron cascade model PACIAE [26[ is 
applied to generate real events. Each real event obeys 
dynamical conservation laws (such as net baryon num- 
ber conservation, energy and momentum conservations, 
etc.), which cause correlations among particles in a single 
event. As the PACIAE model simulation for a nucleus- 
nucleus collision is parallel to the experiment of a nucleus- 
nucleus collision, the NA49 method [13] is also employed 
to generate the mixed events. This means that the mixed 
events are constructed by combining particles randomly 
selected from different real events, while reproducing the 
event multiplicity distribution of real events. We have 
checked that the dynamical conservation laws really do 
not exist in the mixed event and there are only statistical 
fluctuations caused by the effects of finite event number, 
finite size, and the experimental finite detector resolu- 
tions. The "non-statistical moments" of conserved quan- 
tity event distributions are defined to be the difference 
between the moments derived from real events and the 



III. RESULTS 

The net-proton event distribution and the correspond- 
ing higher moments in relativistic Au+Au collisions at 
^Jsmn from 11.5 to 200 GeV are calculated with both 
the PACIAE real events and the corresponding mixed 
events. It is found that the PACIAE results of mean, 
variance, skewness, and kurtosis in Au+Au collisions are 
in agreement with STAR data. 

Shown, as an example, in Fig. [T] are the results of 
PACIAE real events (open symbols) compared with the 
STAR data (solid symbols) of net-proton event distribu- 
tion in 0-5% most central Au+Au collision at ^/sjvjv=39 
GeV. Note that the spectator protons have been excluded 
here. One sees from Fig. [T] that the agreement between 
the STAR data and PACIAE results is satisfactory: with 
two distributions having nearly the same peak location 
and close each other until the half height. The devia- 
tion of the two distributions gets visible around the tails, 
however, the contribution of the tails to the moments is 
expected very small as the corresponding probability is 
rather low. In our knowledge there was not yet compar- 
ison between the STAR data of net-proton event distri- 
bution and model calculations, except the parameter fit 
in H3. 

The centrality dependent STAR data of the moments 
of net-proton event distribution in Au+Au collision at 
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FIG. 4: Energy dependence of non-statistical higher moments 
and their products of the net-proton event distribution in 0- 
5% most central Au+Au collisions. 



Y / s]v7V=39 GeV (cf. first quotation of Ref. [16j|) are com- 
pared with the results of PACIAE real events in Fig. [2] 
In this figure one sees the agreement between the STAR 
data and PACIAE results for M and S. However the 
agreement in a and k is not as good as in M and S. This 
reveals that the tails of the net-proton event distribution 
affect the even moments stronger. 

Fig. [3] shows the energy dependence of the STAR na 2 
data and the results of PACIAE real events in Au+Au 
collisions. The STAR data characteristic feature of al- 
most energy independent na 2 is approximately repro- 
duced by PACIAE within error bars. 

We give in Fig. 2] the PACIAE results of energy de- 
pendence of the non-statistical skewness (panel a), ku- 
torsis (b), Sa (c), and na 2 (d) calculated for the 0-5% 
most central Au+Au collisions in the STAR acceptances 
of \y\ < 0.5 and 0.4 < p T < 0.8 GeV/c according to 
Eqs. ((T][8l). The energies applied in the calculations are 
77^7=11.5, 19.6, 39, 62.4 and 200 GeV, with the to- 
tal number of generated real events being 3.0 xlO 5 , 3.0 
xlO 5 , 2.4 xlO 5 , 1.2xl0 5 and 1.2 xlO 5 , respectively. The 
same number of mixed events are then constructed cor- 
respondingly. It is found in this figure that the non- 
statistical S, k, Sa, and no 1 change sign at ^/snn ~ 60 
GeV. Here a first evidence of singularity is seen in the 
energy dependent S and Sa. These signs, of course, are 
not seen in the results calculated by mixed events. They 
are also not showing up in the results calculated by real 
events and in the STAR data (results of full fluctuations 
indeed) [H, HH , because dynamical fluctuations are al- 
ways submerged in the statistical fluctuations. It seems 
that the higher moment singularity be implicated in the 
non-statistical fluctuation. 
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FIG. 5: (Color online) Energy dependence of non-statistical 
S and Sa of the net-proton event distribution in 0-5% most 
central Au+Au collisions. 



To check the reliability of the sign changes in Fig. 01 we 
have recalculated the energy dependent non-statistical S 
and Sa with 1/2 and 4/5 total number of events and com- 
pared the results with the ones calculated with the full 
number of events, as shown in Fig [5] We see in this figure 
that the convergence is quite good for ^snn= 11.5, 19.6, 
39 and 200 GeV and satisfactory for 62.4 GeV. Thus the 
results shown in Fig. 2] are reliable. The sign changes of 
S and Sa shown in this figure may be a response to the 
prediction in (l2l |. 

The results in Fig. 0J when compared to the estimates 
in the second quotation of [ll| . indicate that the freeze- 
out point is very close to the hadronization point as 
the value of skewness is rather small. Furthermore, at 
lower beam energies the skewness is positive, indicating 
a freeze-out closer to the QGP side in the hadronization 
process. But at higher beam energies the skewness is neg- 
ative, indicating more dominance of the hadronic side of 
the phase transition. The kurtosis becomes significantly 
negative, although small at higher beam energies, indi- 
cating that the freeze-out is indeed in the phase transition 
domain. Since a number of dynamical ingredients are in- 
troduced in the PACIAE (PYTHIA) model, the detailed 
roles and effects of the dynamical ingredients relevant to 
this phenomenon have to be studied later. 



IV. CONCLUSION 

In summary, we have calculated the higher moments of 
net-proton event distributions in relativistic Au+Au col- 
lisions at ^/satat from 11.5 to 200 GeV with real events 
generated by the parton and hadron cascade model PA- 
CIAE [11. The PACIAE results of centrality dependent 
net-proton M, a, S, and n in Au+Au collisions are con- 
sistent with the STAR data. 

We have studied the non-statistical fluctuations of net- 
proton event distributions. The mixed events are con- 
structed according to PACIAE real events and the non- 
statistical moments are calculated as the difference be- 
tween the moment calculated from real events and the 
one from mixed events. The non-statistical S, k, Sa, 
and na 2 appear to change signs at ^Jsnn ~ 60 GeV. The 
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clear sign change in the energy dependent non-statistical 
S and Sa may reflect some singularities. However, this 
has to be confirmed by the STAR Beam Energy Scan 
non-statistical fluctuation data later. 
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